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ABSTRACT 

The performance of routine FNWC sea level pressure anal- 
yses and forecasts in the Beaufort Sea is evaluated for driv- 
ing an ice model. The geostrophic winds determined from the 
FNWC pressure fields are compared to geostrophic winds de- 
termined from the AIDJEX experiment observations. A simu- 
lation with the AIDJEX model has been run using FNWC sea level 
pressure analysis fields for the period 15 to 25 May 1975 
and compared to the AIDJEX results. 

RMS sea level pressure errors of about 3 mb were found 
with analysis errors reduced by about 1 mb when one AIDJEX 
observation was included in the analysis. Geostrophic winds 
had speed differences of about 40% and direction differences 
of about 50°. The FNWC simulation with the AIDJEX model 
produced a position error of about 7 km for one station 
over the ten day run. 

Recommendations are made to improve the FNWC hemispheric 
fields in the Arctic for application to ice forecasting. 
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I. INTRODUCTION 



Recent developments have demonstrated the strategic 
importance of the Arctic to the security of America. Soviet 
Delta-class submarines are on patrol out of their homeport 
of Murmansk armed with nuclear tipped missiles able to strike 
the heartland of America from anywhere within the Arctic Ocean. 
Soviet strategy is to seek open spots in the ice pack as 
launch points for their attacks [Anonymous, 1977]. The 
Northern Sea Route represents a major Soviet strategic inter- 
est [Synhorst, 1973] . As reported in the Monterey Peninsula 
Herald , (1977) , the Soviets have demonstrated their ability 
to penetrate to the North Pole with an atomic ice breaker and 
announced their intention to open a great circle shipping 
route through the Arctic Ocean, events which make their inter- 
est all the more ominous. German submarines in World War II 
demonstrated a threat to the Northern Sea Route by their 
operations in the Kara Sea [Synhorst, 1973] . In the event 
of war any attempts by U. S. submarines to counter the Soviet 
SSBN threat in the Arctic or sever the Soviet lines of com- 
munications via the Northern Sea Route must be aided by re- 
liable knowledge of ice conditions in order to make acoustic 
performance predictions and optimum under-ice track routing. 

The North Slope of Alaska holds a critical supply of 
America's lifeblood, oil. The Navy provides forecasts of ice 
conditions vital to the continued free passage of U. S. 
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maritime resupply vessels to this valuable area. Again, re- 
liable ice forecasts are crucial to America's interests. 

A. NATURE OF THE PROBLEM 

A knowledge of ridge and lead formation and distribution 
are of importance in meeting several Navy requirements. For 
example, the long, smooth refrozen leads are the basic landing 
strips on the pack. Knowing the location and ice thickness 
of these refrozen leads would be of great help in carrying 
out scientific research and logistic support on the pack. In 
shallow areas such as the Chukchi Sea the ice keels beneath 
ridges can form a considerable hazard to under-ice navigation. 
Ice free leads or thin ice can act as places for emergency 
surfacing of a submarine. Penetration of the ice is impeded 
when the pack is under high internal stress. Combined with 
a knowledge of ice thickness, a knowledge of the stress state 
of the ice would aid in these decisions, and weapons and 
sensor deployment through the pack. A submariner might find 
accurate forecasts of these conditions very reassuring on a 
submerged transit of the Arctic. 

The two features in the polar environment that most 
strongly influence underwater sound are the permanent ice 
cover and the sound speed structure of the water. The essen- 
tially positive sound speed gradient causes sound rays to be 
generally refracted upward and then reflected and/or scattered 
from the rough ice-water interface [Diachok, 1976] . Both 
high and low frequencies are rapidly attenuated, the former 
by scattering losses from the ice cover, the latter by the 



10 



fact that very low frequencies are not effectively trapped 
in the channel [Urick, 1975]. Recently it has been demon- 
strated that back-scattering strength is controlled virtually 
exclusively by the size and number of sea-ice ridges [Diachok, 
1976] . 

The principal source of ambient noise in the Arctic and 
sub-Arctic water is large scale deformation of ice floes 
[Diachok, 1976]. When the ice is not continuous, as near the 
ice edge, noise levels 5 to 10 dB higher than those measured 
at the same sea state in ice-free water have been observed 
[Urick, 1975]. The mechanisms responsible are related to 
wind-generated ocean wave interactions with sea ice at the 
boundary [Diachok, 1976]. Under a continuous ice sheet with 
calm winds and rising temperatures very low levels of ambient 
noise, well below the lowest Knudsen curve for sea-state zero, 
have been found. On the other hand, levels some 40 dB higher 
have been observed when the ice had been cracking under fall- 
ing air temperature [Urick, 1975] . The acoustic propagation 
and ambient noise factors are both important in under-ice 
navigation, submarine tracking and localization, and weapons 
performance. 

In addition to its effect on the under-ice environment 
the ice cover has a significant effect on the climate of the 
Arctic by acting as a lid on the ocean. A typical mesoscale 
region in the ice pack can be expected to contain open water, 
young ice tens of centimeters in thickness, perennial and 
deformed ice a few meters thick. Of concern in atmospheric 
modeling is the fact that open water and thin ice can exert 
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a large influence on regional heat exchange rates. For 
example , results by Maykut [1976] indicate ice 0-40 cm thick 
allows a net heat exchange to the atmosphere one to two orders 
of magnitude greater than perennial ice a few meters thick. 

Ice thickness greater than 100 cm exhibit little change in 
heat flux as the ice thickens. In the Central Arctic during 
winter , estimates of open water and young thin ice are less 
than 1% and 8 to 12% , respectively. Thus , even at this small 
percentage the effects of these two ice conditions can begin 
to dominate the large scale heat exchange and alter boundary 
conditions. Atmospheric models which assume the ice cover 
is a continuous three meter layer may seriously underestimate 
the heat input to the boundary layer. 

At the present time the Navy ice prediction models are 
only intended to serve as an aid to hand analysis. The wind 
drift model used by the Navy ice forecasting group at the 
Fleet Weather Facility, Suitland, Maryland, gives ice speed 
as a linear function of surface wind speed [Skiles, 1968], 
Moving to the right of the wind, the drift angle is given as 
an inverse exponential function of wind speed. Under de- 
velopment is a model to provide daily wind drift after 
Shuleikin [1953] and growth after Zubov [1945] at 62 coastal 
weather reporting stations [Gerson and Simpson, 1976, and 
Gerson, 1975], Barnett [1976] has presented a scheme for 
estimating the general ice conditions during the shipping 
season of the North Slope sea route using the height of the 
1000 mb surface at two geographic locations (52N, 100E, and 
7 ON , 140E) . 
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All these models are conceptually extremely simple in 
the light of our present understanding of the complexity of 
sea-ice dynamics and they will probably not be capable of 
meeting the expanding needs of the Navy for arctic ice fore- 
casting. Several more complete dynamical models of the arctic 
ice pack have been formulated [Campbell, 1965; Coon, et al , 
1977; Hibler, 1977]. Therefore, a new model is in order, 
but first an assessment of present Navy analysis and forecast 
fields in the Arctic is needed to determine if these are 
capable of meeting the needs of the more complex models. 

B. DESCRIPTION OF THE ARCTIC ICE PACK 

The Arctic ice pack forms a lid over most of the North 
Polar Sea, constantly changing its extent of coverage, 
patterns of motion and thickness distribution. The seasonal 
range of ice coverage is illustrated in Figure 1. Of note 
is the extreme variability of the position of the ice edge, 
resulting in a 20% change in the areal coverage, receding in 
summer far enough to allow a shipping season along the north 
coast of Siberia and the North Slope of Alaska. 

The summer season brings melting to the central pack, 
reducing the average thickness by about one meter and pro- 
ducing large open water areas [Orvig, 1970] . These open 
water areas effectively reduce the large scale strength of 
the ice pack to zero and allow the floes to move freely. 

While the day-to-day movement of individual ice floes 
resembles a random walk, the long term general circulation of 
the ice pack has a regular pattern. Two dominant features 
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are a trans-polar drift from the East Siberian Sea, across 
the North Pole, and out of the Arctic into the East Greenland 
Current, and an anti-cyclonic gyre in the Beaufort Sea. These 
large scale climatic motions resemble the atmospheric climatic 
circulation patterns in the same way that the lower latitude 
ocean circulations resemble climatic wind patterns [Campbell, 
1965] . 

In contrast to the relatively free movement of the ice 
in the central pack, the coastal ice remains attached to 
land in winter, preventing longshore movement. The transition 
from zero velocity of the land fast ice to a velocity of a 
few centimeters/sec in the central pack occurs in the shear 
zone, normally found a few kilometers offshore where sharpe 
ice velocity gradients occur. 

The central pack has an average thickness of about three 
meters and it is made up of individual floes with dimensions 
of a few meters to several hundred meters [Orvig, 1970]. In 
winter these floes are separated by ridges up to several 
meters high and interlaced with leads of open water up to 
0.5 km wide and several kilometers long, leads which quickly 
refreeze to form large smooth areas of young ice which even- 
tually form ridges and new floes. These ridges and leads are 
formed by the relative motion of the ice whose principal 
driving force is the wind. The nature of ridging is extremely 
complex. Reduced to a simple description, sea-ice ridges are 
rubble piles formed by collisions and shear interactions 
between adjacent floes. The ridges above level ice are 
accompanied by keels below the surface, keels which may 
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extend many times as deep as the ridge is high. Since on a 
large scale the ice is densely fractured it has little tensile 
strength and readily opens when subjected to a stretching or 
divergent strain [Pritchard and Colony, 1976] . 
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II. HISTORICAL BACKGROUND 



The modeling of ice dynamics can be divided into three 
categories; empirical modeling, small-scale analytical model- 
ing and synoptic-scale analytic modeling. The empirical 
modeling began with the observations by Nansen during the 
drift of the Fram across the Arctic Ocean (1893-96) , who 
observed that ice drift speed was about 2% of the wind speed 
and generally to the right of the wind direction at an average 
angle of 28°. Sverdrup, drifting with the ice pack in the 
Maud , (1919-25), found a wind factor for ice of 0.014 to 
0.024 and a drift angle from 18° to 40°, the latter when the 
ice was loose in summer. 

Zubov [1945] further developed the empirical relationship 
of wind and ice drift. His analysis of the drift of the 
Sedov yielded a wind factor of 0.015 and drift angle of 29°, 
which compare favorably with Nansen. Zubov further postulated 
that ice drift follows the sea-level isobars, i.e., it is 
parallel to the sea level geostrophic wind because the counter- 
clockwise turning of the wind through the atmospheric friction 
layer is approximately equivalent to the ice drift angle to 
the right of the surface wind. 

Early work in analytic modeling of small-scale ice 
dynamics was done by Shuleikin [1953] who derived the force 
balance and motion of the ice due to the wind stress. However, 
his theories were limited to 10/10 ice coverage. Fukutomi 
[1948] treated ice drift by considering the effects of ice 
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concentration, roughness and ocean currents in an analytical 
model . 

Modern high speed computers have made it practical to 
develop synoptic-scale ice drift models which span a broad 
range in complexity. Simplest are the wind drift models 
using empirical relationships. An example of the use of wind 
drift relationships is the numerical model developed by Knodle 
[1964], which combines Shuleikin's and Fukutomi ' s theories, 
yielding a wind drift field which is a function of ice con- 
centration for Baffin Bay. The most complex are analytical 
models which consider the interactions within the pack ice. 
Because of the internal stresses transmitted through the ice 
pack a solution which neglects the interaction of the ice 
cannot completely describe the dynamics. 

The first models which attempted to include the inter- 
action of the ice assumed it to have a viscous nature. In 
formulating his ice circulation model Campbell [1965] considers 
individual ice floes as "... fluid elements with an average 
area which is small compared with the ocean area, [thus] the 
ice sheet may be viewed as a film of highly viscous fluid 
suspended between two less viscous fluids, air and water." 

While Campbell and several other authors have produced viscous 
models which describe the general circulation of ice in the 
Arctic, all have limitations, the most notable is their in- 
ability to model the coastal shear zone. 

Hibler [1977] has developed a viscous model which simu- 
lates large scale (100 km) sea-ice dynamics; the model has 
been run for simulations 8 years in length. Since this model 
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is not hampered by numerical stability problems it is able 
to utilize large time steps (1 day) and thus perform long- 
term simulations, while still retaining the ability to be run 
in smaller time and space scales. This model incorporates 
an ice thickness evaluation which combines mean thickness 
and compactness. Despite his model's success Hibler states, 
"... particular elastic-plastic rheologies may well be more 
useful in reproducing many detailed small scale effects." 

One of the most recent models, which takes a departure 
from viewing the ice as viscous in nature, is that formulated 
by the Arctic Ice Dynamics Joint Experiment (AIDJEX) program. 
It models the ice as an elastic-plastic material whose yield 
strength depends on the amount of thin ice present. A major 
difference between the plastic and fluid models is that the 
plastic material responds independently of the rate of the 
deformation [Coon, Hall and Pritchard, 1977]. In simulations 
covering periods of about 10 days the AIDJEX model has pro- 
duced good results in describing both large and small scale 
features of ice dynamics and depicting the coastal shear 
zone [Coon, et al, 1977]. Because of the potential success 
of the AIDJEX model in addressing these important Navy re- 
quirements for ice forecasting, it was selected as a means 
of evaluating the Navy's atmospheric analysis and prognostic 
fields for use in an ice dynamics model. 
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III. DESCRIPTION OF THE AIDJEX MODEL 



The AIDJEX ice model consists of three major parts: the 

momentum equation, a set of constitutive laws which describe 

the stress-strain behavior of the ice in the horizontal plane, 

and the ice thickness distribution [Coon, et al, 1976]. 

The first element of the model is the momentum equation 

which states that the acceleration of an element of ice of 

mass (m) is related to the air stress, (r ) , the water stress 

3. 

(r w ) , the coriolis force, and the divergence of internal 
stresses ( 7 . a ) . The momentum equation is written: 



V • a 



■mf kxv 
c 



mf kxv 
c g 



( 1 ) 



where 



v 



-mf kxv 
c 

f 

c 

mf kxv 
c g 

v 

g 



= ice velocity 
= coriolis force 
= coriolis parameter 
= sea surface tilt force 
= geostrophic current 



Air stress is computed as a quadratic function of the 
geostrophic wind [Pritchard, et al, 1976]: 



? a “ C D I 2 ' ?a 5 (2) 

where 

p = air density 

cl 

C D = geostrophic drag coefficient 
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u 



B 

- a 
a 



= geostrophic wind velocity 
cosa -sin a 

sina cos a 

= counterclockwise drag angle 



Water stress is represented by a quadratic drag law 

similar to that used in the atmospheric boundary layer. The 

stress on the lower surface of ice (r ) is a function of the 

w 

ice velocity (v) relative to v : 

g 



where 




v g | B (V 



V 



(3) 



P w = water density 

C w = water drag coefficient = 5.5 x 10 



-3 



B = 



cos (7r-/?) -sin (7 r-/?) 

sin(7r-/?) cos (7 r-/?) 
yS = drag angle = 23° 

Units are in the CGS system. The geostrophic flow is taken 
to be given by the long term observed values [Pritchard, 
et al, 1976] . 

The stress deformation constitutive law has been chosen 
as an elastic-plastic response. A plastic model possesses 
many of the properties thought to be desirable in describing 
the response of ice to wind stress. Specifically, it is able 
to simulate rapid variations which approximate discontinuous 
behavior in nearshore regions and simulate flow resembling a 
viscous material when a large enough load is applied. The 
constitutive law relates internal ice stress ( a ) to the 



20 



large scale strain in the ice sheet. The ice is a stiff 
elastic material allowing elastic deformation of only about 

0. 1% before reaching a yield strength in compression and 
shear. Stresses are restricted to have no positive principal 
values because the material is densely fractured and will not 
support stress in any opening mode [Colony, 1976] . 

The thickness distribution measures the percentage of 
area covered by ice in each differential range of thickness 
within each element of the grid. A typical example is shown 
by the solid line in Figure 2. The thickness distribution 
determines the mass per unit area of the element as well as 
its mechanical properties. The thickness distribution is 
affected by both thermodynamic and mechanical processes. 

For example, the melting of a few centimeters off the top 
surface would shift the thickness distribution to the left, 

1. e., to the dashed line in Figure 2. During the winter, all 
but the very thickest ice grows thicker, shifting the curve 
to the right. When the element converges slightly, the thin 
ice is ridged into thick ice, causing the shape of the distri- 
bution to change as some thin ice disappears and a smaller 
increment of new thick ice is produced, as given by the 
dotted line in Figure 2. If the element diverges, open water 
forming between the floes would be shown in the distribution 
by a delta function at zero thickness. As time passes this 
open water would freeze into thin ice, moving the delta function 
to the right [Coon, et al, 1976] . In the AIDJEX model the 
growth rate is defined as the monthly climatological average 
[Pritchard and Colony, 1976]. 
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One of the desirable properties of the AIDJEX model is 
the potential to express yield strength as a function of the 
thickness distribution. However, thus far, a satisfactory 
relationship has not been found and the yield strength is 
varied as a unique input parameter. The redistribution 
function and the energetics argument that enables strength 
to be determined from the thickness distribution will be 
reformulated as further work provides direction [Pritchard, 
et al, 1977] . 

The computational grid with manned camps and buoy loca- 
tions for 15 May 1975 is shown in Figure 3. The grid was 
assembled according to the following rules: (1) The boundary 

may be a coastline or a dynamic boundary determined by the 
buoy positions, etc. Positions of the data buoys and the 
manned drifting stations were found using the Navy Navigational 
Satellite System and the buoy ring motion determines the 
boundary velocity conditions. This is the only boundary con- 
dition used in the model. (2) All interior cells must be 
quadrilaterals and all cells next to the boundary must be 
quadrilaterals or triangles. To minimize computer time, 
linear dimensions of all cells were made approximately equal. 
The stability of the difference scheme depends on the Courant 
restriction, C \ , which means that the time step is pro- 
portional to the cell dimension, specifically the dimension 
of the smallest cell. The factor, C, is a function of ice 
strength; as ice strength increases the time step must de- 
crease . 
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Initial conditions for a model run are contained on a 
computer file, magnetic tape or disk, which must be generated 
from a history of ice motion derived from buoy and manned 
camp positions. This file, called a DUMP file, contains 
positions of the grid points, initial velocities of grid cells, 
stress state, and ice thickness. Periodically during the 
model run this file is regenerated, providing a set of initial 
conditions to restart the run. Thus the model can be stopped 
and restarted several times during a lengthy run to allow 
checking of the solution or to rerun the problem and look at 
some event on a finer time scale. For example, if the solution 
were being saved for plotting every four hours, it might be 
desirable to return and look at the solution every twenty 
minutes for some four-hour period. The DUMP file also serves 
as the data source for plotting two-dimensional fields. 

One additional output from the model is a time history 
of data from 20 selected nodes. Saving the solution at only 
a fraction of the nodes enables one to plot the solution 
densely in time, enabling high frequencies to be seen without 
using excessive amounts of storage [Don Thomas, AIDJEX, 
personal communication] . 

The parts of the model, i.e., the momentum equation and 
the constitutive laws are solved simultaneously. Detailed 
equations and methods of solution are given by Pritchard 
and Colony [1976] . 
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IV. EVALUATION OF FNWC SEA-LEVEL PRESSURE FIELD AND 
GEOSTROPHIC WINDS 



Fleet Numerical Weather Central (FNWC) , located in 
Monterey, California, is responsible for providing all numer- 
ical environmental forecasts in support of Naval operations. 
FNWC is linked to the National Meteorological Center (NMC) , 
the Air Force Automated Weather Network (AWN) , and in addition 
receives observations from Navy units throughout the world 
[U. S. Naval Weather Service, 1975]. The numerical products 
produced at FNWC include sea-level pressure, upper-air analyses 
and prognoses, marine layer winds, sea-surface temperature, 
and many others produced on a polar stereographic projection 
from data analyzed on a rectangular grid. Most fields are 
resolved on a 63 x 63 point Northern Hemispheric grid but 
other grids are used for some products. Atmospheric prognoses 
(forecasts) are made twice a day at 0000GMT and 1200GMT for 
each 12-hour period out to 72 hours. Many computational 
algorithms are used to generate the numerical products but 
one of the most important is that which blends available 
information to formulate the best possible sea-level pressure 
field. This is known as the Fields by Information Blending 
(FIB) technique. 

A. DESCRIPTION OF THE FIELDS BY INFORMATION BLENDING (FIB) 

TECHNIQUE 

At FNWC, sea-level pressure is analyzed to a 125 x 125 
polar stereographic projection hemispheric grid (grid spacing 
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190 km at 60 °N) using the FIB technique [Holl and Mendenhall, 
1972]. While analyzed on a 125 x 125 grid the resultant field 
is archived on a 63 x 63 grid which has a 381-km spacing at 
60 °N . 

The FIB generation of the sea-level pressure analysis is 
illustrated in Figure 4. The basic steps in the FIB analysis 
technique are : 

1. The 6-hour old analysis is extrapolated to the 
present analysis time, assembled at each grid point 
and combined with the 6- or 12-hour old prognosis 
to produce a first guess field. 

2. New wind and pressure reports are assembled, re- 
jecting errors and weighting individual values 
according to their reliability. 

3. Available information fields are blended into a 
modified first guess field. 

4. Fields of sea-level pressure, sea-level pressure 
gradient and the Laplacian of the pressure field 
are analyzed and blended to obtain the best fit 
field. 

5. New reliability weights are determined. 

6. A new gross error check and reevaluation of data 
is made. 

7. Assembly, blending and resultant weight operations 
are recycled. 

8. The final best-fit pressure field is output. 

Since the resulting pressure field is a blend of avail- 
able information, in the absence of observations the analysis 
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will be largely based on the previous prognosis. However, 
the quality of the prognosis is dependent upon the accuracy 
of the previous analysis. Thus, in data sparse areas such 
as the Arctic, the analysis is prone to greater error than 
a data rich area. 

B. SEA- LEVEL PRESSURE EVALUATION METHOD 

Evaluation of FNWC 1 s atmospheric sea-level pressure 
analysis was performed by using AIDJEX sea-level pressure 
observations as a standard. Twenty-day files of sea-level 
pressure observations and geostrophic wind analyses were 
obtained from AIDJEX on magnetic tape. The period from 1 May 
through 19 July 1975 (AIDJEX days 121-200, four 20-day inter- 
vals) was selected for detailed evaluation.'*' AIDJEX data 
represent high density, high quality observations from several 
manned camps and remote buoy stations whose positions on 
18 May 1975 are shown in Figure 5. The FNWC sea-level pres- 
sure analysis for this time is shown in Figure 6. 

FNWC ' s 6-hour ly analyses and 12-hourly prognoses of sea- 
level pressure fields were interpolated to each AIDJEX ob- 
servation station. The RMS pressure difference (where pres- 
sure difference is A p=p pnwc _P AIDJEX^ at AIDJEX stations 

in each 20-day interval was calculated for the analysis, 24-, 
48-, and 72-hour prognoses, and the results are summarized 
in Figure 7. The RMS pressure difference (a measure of the 



^AIDJEX days are Julian days measured from 1 Jan 1975 
and continue to be consecutively numbered into 1976. Thus 
2 Jan 1976 is AIDJEX day 366. 
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quality of the product) increases from 2.3 mb in the first 
interval iA, 1-20 May) to 3.4 mb in the last interval (D, 

30 June to 19 July) . For the prognostic times the error 
increases by about 0.6 mb for the first 24-hour period and 
thereafter by about 2 mb for each 24-hour projection. For 
some reason the second and third time intervals (B , 21 May to 
9 June and C, 10-29 June) show an increase in quality over 
the first interval. However, the quality is always poorest 
in the final interval which is the summer period. 

An explanation for the larger RMS A P values in the 
fourth time interval is found in Figure 8 which shows a time 
series of atmospheric sea-level pressure at the AIDJEX main 
camp over the entire 80-day period. The greatest variability 
of the pressure is in the fourth period implying an increased 
synoptic activity level. This increased activity combined 
with lack of observations probably accounts for the decreased 
quality. 

An attempt was made to determine if and when observations 
from AIDJEX were included in the FNWC analysis. A search was 
made of FNWC archives for any reports submitted by AIDJEX. 

A total of 95 AIDJEX reports were found in the archives, re- 
ported only from the main AIDJEX camp. However, due to the 
nature of the archiving system late reports (up to 24 hours) 
are archived and it is not possible to determine if any par- 
ticular report was included in the analysis. RMS A P for 
analyses corresponding to archived reports is shown by the 
line marked R on Figure 7. A careful examination of Figure 
8 reveals the effect of the inclusion of an AIDJEX pressure 
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observation in the FNWC analysis. Time of archived AIDJEX 
observations are indicated by "tick" marks near the lower 
edge of the graph. The circled areas illustrate a few times 
corresponding to archived observations when the analysis has 
made a sharp change in trend to become coincident with the 
value of the observations. It can be inferred that the AIDJEX 
reports did reduce errors in the FNWC analysis by 0.7 mb on 
the average and that any observations from the central Arctic 
will probably result in an improved analysis. 

Walsh [1977] has calculated the RMS height errors for 
the 1000 mb height forecasts of the NMC primitive equation 
model at 16 grid points covering roughly the same area as 
the AIDJEX observation stations. His technique involved 
finding the height error of the 24-, 48-, and 72-hour fore- 
casts using the verifying analysis as a standard. When con- 
verted to millibars (8 m = 1 mb) these errors range from 
2-3 mb at 24 hours to 5-6 mb at 72 hours. This technique of 
using the verifying analysis vice observations as a standard 
was applied to the FNWC data, again interpolated to the same 
AIDJEX station locations. The RMS A P increases from 2 . 5 mb 
at 24 hours to 6.5 mb at 72 hours. These are comparable to 
the results found by Walsh. However, they show a marked 
improvement over the use of observations as a standard, the 
improvement ranging from 1.0 at 24 hours to .3 mb at 72 hours. 
This suggests that the use of the verifying analysis as a 
standard for forecast evaluation yields optimistic results 
primarily because the analysis is largely based on the pre- 
vious prognosis in the Arctic area. 
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Sufficient pressure data were available from AIDJEX to 
evaluate FNWC analyses from 1 May 1975 through 24 April 1976. 
The results of the evaluation are shown in Figure 9. A P 
began at about 0.3 mb (higher pressure than actually existed) 
and continued increasing into the summer until the end of 
August when the value reached 2.7 mb. It decreased as fall 
approached becoming negative in mid-October, and generally 
remained negative throughout the winter, becoming positive 
the following March. 

The RMS A P remained high throughout the summer, de- 
creased slightly in the fall and remained lower than the sum- 
mer values throughout the winter. The sharp increase in 
March can be attributed to increased synoptic activity. An 
evaluation of whether these fluctuations of RMS A P and A P 
follow an annual cycle will become possible when a longer 
record of complete observations in the Arctic become avail- 
able . 



C. METHOD OF EVALUATION OF GEOSTROPHIC WINDS 

An evaluation of FNWC geostrophic winds was performed to 
determine their suitability for air stress computation. In 
the AIDJEX model and probably any other operational model the 
wind stress is determined from the geostrophic wind. There- 
fore, the quality of the sea-level pressure field is an im- 
portant consideration. The AIDJEX sea-level pressure field 
is expressed by a least squares polynomial fit to all data 
points [AIDJEX, 1976]. However, geostrophic winds are de- 
pendent on the pressure gradient rather than the actual 
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pressure values at any point. Geostrophic winds, as obtained 
from AIDJEX, are found by differentiating the polynomial ex- 
pression along grid directions at each observation station. 

To obtain FNWC geostrophic winds, a fourth order finite 
differencing technique was applied to the FNWC grid point 
sea-level pressure value to determine the pressure gradient 
components. FNWC geostrophic wind components were then inter- 
polated to AIDJEX observation stations and compared with 
AIDJEX geostrophic winds. The RMS wind speed and direction 
differences were calculated and the results are summarized 
in Figure 10 where wind speed differences have been normalized 
by the AIDJEX value. Aggregate results of all four evaluation 
periods are marked by *G which shows increased differences 
as the forecast period increases. In particular the direction 
differences increase from about 50° to 75° from analysis to 
the 72-hour forecast, respectively. The fractional speed 
errors show an increase from 0.65 to 1.35 for the same periods. 

If only wind speeds in excess of 5 m/sec are considered, 
the RMS differences decrease significantly, and only a slight 
increase from 0.4 to 0.5 is found as the forecast period 
lengthens. These values are of particular interest since the 
air stress is taken as proportional to the square of wind 
speed and the higher wind speeds will be of greater importance 
in determining ice drift. There is a relative minimum in both 
direction and speed differences in the 21 May to 9 June time 
period which corresponds to the reduced RMS pressure error 
in this time frame. However, for some reason the direction 
and speed difference for period C do not similarly correspond 
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with the relatively low RMS pressure errors for that time 
frame. 

D. SIMULATION RESULTS 

To more completely evaluate the usefulness of the FNWC 
sea-level pressure field in driving an ice dynamics model, 
a simulation using the AIDJEX model was made. The AIDJEX 
model was programmed on the FNWC 6500 computer to simulate 
ice movement during the period 0300GMT 15 May to 1200GMT 
25 May 1975. Air stress derived from the FNWC sea-level 
pressure analysis field was used in the simulation while all 
other input parameters and boundary conditions were specified 
as those used by AIDJEX. The air stress field was computed 
using equation (2) and assembled in a form compatible with 
the AIDJEX model using the following values of air density 
( P ) , drag coefficient (C n ) and turning angle ( a ) : 

3 D 

= 0.00135 gm/cm 3 
C D = 11 x 10" 4 
a = 20 ° 

C D and a are the same values used in AIDJEX model simulations 
[Eric Levitt, AIDJEX, personal communication] . 

Figures 11 and 12 show time series plots of both FNWC 
and AIDJEX computer air stress magnitude and direction at two 
grid points. It is clear that the FNWC air stress is smaller 
than the AIDJEX air stress, in particular for the important 
time periods days 135-136 and 137—140, when there is a high 
air stress. At node 7,7, which corresponds to the main 
AIDJEX camp, during day 135 the FNWC air stress is 20% of 
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AIDJEX values. From day 137 to 140 the FNWC air stress is 
in general much less than AIDJEX. The values at node 3,7 
follow a similar pattern, the FNWC stress being a small 
fraction of the AIDJEX value for the first day, and the FNWC 
stress generally retaining the trends of the AIDJEX stress, 
but at reduced levels. In addition to the magnitude differ- 
ences there are direction differences of up to 150°, however, 
these are generally limited to the light wind conditions. 

The direction differences for the critically high stress 
period 137-140 are small, i.e., generally less than 30°. 

Figures 13 and 14 show vector field plots of air stress 
at 1200GMT 18 May for the AIDJEX and FNWC model runs, respec- 
tively. The AIDJEX plot shows considerable spatial variation, 
with an area of high stress over the southern portion of the 
region. The FNWC stress field plot shows less spatial vari- 
ation and significantly smaller air stress over the north- 
eastern portion of the region. 

The trajectories computed using FNWC and AIDJEX stress 
fields for the main AIDJEX camp are shown in Figure 15. The 
total displacements of the FNWC projected movement is compar- 
able to the AIDJEX movement with only a 3 km difference in 
the final positions. In addition to the comparable total 
displacement, the small-scale features of the AIDJEX movement 
are retained in the FNWC movement. Figure 15 also includes 
the actual movement of the AIDJEX main camp as calculated 
from AIDJEX measured positions [Thorndike and Cheung, 1977] . 

It is clear that neither AIDJEX nor FNWC followed the movement 
exactly. During the first several days the FNWC movement 
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more closely resembles the true motion, while in the final 
positions the AIDJEX simulation is closer, namely 3 km versus 
4.5 km. 

A more detailed analysis is possible by looking at the 

x- and y-position time series for the main camp. Figure 16 

shows the x- and y-position time series relative to the AIDJEX 

coordinate system [Thorndike and Chueng, 1977]. It is evident 

that the error in the y-position was accumulated largely in 

the first two days. Referring back to Figure 11, the FNWC 

2 

air stress is only about 0.3 dynes/cm compared to 1.5 dynes/ 
cm^ for AIDJEX. This low FNWC value is the cause of most of 
the difference in the final position. 

The relative motion of two or more stations is important 
for those processes such as deformation which depend on 
differential motion. The trajectory of node 3,7, located in 
the left hand portion of the grid, is shown in Figure 17. It 
has a larger scale of motion than the main camp and the final 
error is in the opposite direction, being about 3 km in the 
negative y direction. This quantitatively implies a larger 
deformation difference between the two stations than the 
absolute motion difference of either single station. Node 
3,7 had no position measurements to allow checking of the 
results of the simulations. 

Figure 18 shows the x- and y-position time series for 
node 3,7. Unlike node 7,7, whose error was accumulated 
largely in the first day, this node accumulated its error in 
the 138-139 day time frame. For node 3,7 the error steadily 
increases from day 137-139.5, this time period corresponding 
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to the maximum difference in air stress between the two runs 
(Figure 11) . 

Figure 19 is a vector plot of the ice velocity field 
for the AIDJEX simulation at 1200GMT 18 May. Recalling that 
the boundary velocity is specified, it is evident that the 
discontinuities in the velocity field along the southwest 
and northeast boundaries represent errors in the calculated 
velocities. Figure 20 is a vector plot of the stress diver- 
gence field at the same time showing relatively large stress 
divergence forces along the south and northeast boundaries. 
The turning of the vector along the south boundary implies 
that the stress is an artifact of the discontinuous velocity 
field and probably has the effect of directing the field 
velocity into agreement with the boundary velocity. 

Figure 21 is a vector plot of the velocity field for 
the FNWC simulation at 1200GMT 18 May. The velocity dis- 
continuity along the northeast boundary is similar to AIDJEX 
but there is no apparent discontinuity along the south bound- 
ary. The regular pattern of the stress divergence shown in 
Figure 22 appears to be a more natural occurrence than does 
the AIDJEX stress divergence field, implying that at least 
for this one time the FNWC simulation is a better representa- 
tion than the AIDJEX simulation of the actual ice dynamics. 

The relationship of the two computed velocities is shown 
in Figure 23. The only large differences in the velocities 
are for days 135 and 137. These correspond to the large 
differences in air stress values. The remaining time period 
shows a high correlation between the two data sets, an even 
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higher correlation than would seem to be supported by the 
significant differences in the air stress alone. The reasons 
for this high correlation can be attributed to the influences 
of the boundary velocity condition. The ice is largely 
affected by internal ice stress during the light wind con- 
ditions and is primarily driven by the boundary velocity 
inputs. This is amply illustrated in Figures 24 through 27 
where the FNWC and AIDJEX air stress fields are radically 
different yet the velocity fields are nearly identical. The 
model is being run by the boundary at this time. 
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V. CONCLUSIONS 



A comparison of FNWC sea-level pressure analyses and 
prognoses with observations made by AIDJEX during the spring 
and summer of 1975 shows: 

1. A RMS error of approximately 3 mb for the analysis 
time increasing to approximately 7 mb at 72 hours. 

2. The RMS error is highest during the summer months, 
increasing from 2.3 mb in May to 3.8 mb in July. 

Comparing the FNWC sea-level pressure analyses and AIDJEX 
observations for the 360-day period from 1 May 1975 to 24 
April 1976 shows: 

1. The RMS error remains 'relatively high throughout 
the summer, decreases into fall and winter with a 
sharp rise in February. 

2. The average 4? increases to a maximum of 1.7 mb in 
August, becoming negative in October reaching a 
minimum of -1.0 mb in November, and then becomes 
positive again in late winter. 

Of particular note is the fact that the inclusion of just 
one observation for the AIDJEX ocean improves the pressure 
analysis by nearly 1 mb. 

A comparison of FNWC generated geostrophic winds with 
AIDJEX generated geostrophic wind reveals: 

1. The normalized RMS speed difference for all wind 
speeds increases from 0.65 to 1.35 from analysis 
to the 72-hour prognosis. 
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2. Considering only wind speeds greater than 5 m/sec, 
the normalized RMS speed difference is lower and 
increases from 0.4 to 0.5 from analysis to 72-hour 
prognosis . 

The movement of two AIDJEX grid nodes using FNWC sea- 
level pressure fields to derive wind stress in a ten-day simu- 
lation closely approximates the movement calculated by AIDJEX. 
Final position errors compared to measured positions of the 
main camp were 4.5 km for the FNWC simulation versus 3 km for 
the AIDJEX simulation. The other node had a difference of 
about 3 km between AIDJEX and FNWC simulations. The majority 
if not all the final position error was accumulated during 
periods of high air stress when FNWC air stress was 20% to 
80% of the AIDJEX value. 

During periods of low wind stress the AIDJEX model is 
largely driven by the boundary velocity conditions as evidenced 
by the high correlation of AIDJEX and FNWC calculated veloci- 
ties during these periods. 

The performance of the AIDJEX model is highly dependent 
on the accurate specification of the boundary velocities. 

Any attempt to use the model as a predictive tool will be 
limited by the use of the boundary conditions. The present 
method is not realistic for forecasting and must be replaced 
by some form of wind driven boundary. A significant degrada- 
tion of the performance of the model can be expected to result 
from the use of a wind driven boundary. 
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Another limitation of the model's use for forecasting 
is the amounts of computer time required for a simulation. 

A 72-hour simulation required about 1.5 hours computer time 
on a CDC 6500 computer. This is fully one-fourth the time 
required to generate a 72-hour 7-level primitive equation 
meteorological prognosis for the entire Northern Hemisphere. 

The important ice thickness distribution was not investi- 
gated in the simulation because it was deleted from the model 
for reasons relating to its contribution to the ice strength 
calculations. However, it may be included later if one 
desires it for other than ice strength considerations. 
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VI . RECOMMENDATIONS 



It is recommended that the Navy continue work on opera- 
tional ice dynamics models. Specific areas of investigation 
should include a formulation of boundary conditions which 
are compatible with a forecasting vice hindcasting scheme, 
preparation of a model which works well in the summer ship- 
ping season, and incorporation and evaluation of an ice 
thickness distribution technique which is independent of ice 
strength calculations. 

Furthermore, every effort should be made to provide 
synoptic sea-level pressure and air temperature observations 
from the central Arctic. The satellite telemetry buoys 
developed and used in conjunction with the AIDJSX field 
experiment offer one of the simplest and most direct means 
of improving the environmental forecasts in the Arctic. With- 
out the real time availability of these observations, good 
predictions of ice conditions encompassing all Navy require- 
ments will remain an impossibility. 
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APPENDIX A FIGURES 
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Extent of sea ice in the Northern Hemisphere, 
and Van Allen, 1971] 
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A representative thickness distribution undergoing changes 
in the text. [Coon, et al, 1976] 
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Figure 3. The AIDJEX computational grid at 0300GMT 15 May 1975. 
Observation stations are marked as : remote buoys ( # ) , manned 

camps ( ■ ), main AIDJEX camp (node 7,7) ( ■ BB) . Node 5,7 is 
marked as ( A ) • Irregular configuration of the grid is a result 
of applying the rules for generating the grid as described in the 
text . 
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Figure 4. Simplified flow diagram of FNWC’s Fields by Information 
Blending (FIB) Sea- Level Pressure analysis technique. 
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Figure 5. AIDJEX sea-level pressure analysis for 1200GMT 18 May 1975. 
Values are sea- level pressure minus 1000 mb. AIDJEX observation 
stations are marked by ( ■ ) . Station marked B3 is the location of 
AIDJEX main camp (Big Bear) . 
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Figure 6. FNWC sea- level pressure analysis for 1200GMT 18 May 1975. 
63x63 grid points are marked ( # ) . Station marked ( £2 BB) is location 
of AIDJEX main camp. Values are sea- level pressure minus 1000 mb. 
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Figure 7. RMS pressure error for FNWC sea- level pressure analysis 
and prognoses at all AIDJEX stations. Analysis time is T=0, T=24 is 
24 hour prognosis, T=48 is 48 hour prognosis, T=72 is 72 hour prog- 
nosis. Dates are shown by A, 1 May to 20 May; B, 21 May to 9 June; 

C, 10 June to 29 June; and D, 30 June to 19 July. Points marked *G 
are aggregate errors for all four time periods. Line marked R is RMS 
error for only those analysis times corresponding to FNWC-archived 
AIDJEX observations. 
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Figure 8a. Time series of AIDJEX sea-level pressure observations 
(solid line) and FNWC analysis (dashed line) for AIDJEX main camp for 
period 1 May to 20 May 1975. "Tick" marks indicate analysis times 
corresponding to FNWC-archived AIDJEX observations. 
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Figure 8b. Time series from 21 May to 9 June 1975. Circled portions 
indicate FNWC- archived AIDJEX observations which obviously had the 
effect of bringing the FNWC analysis to the AIDJEX observed values. 
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Figure 8d. Pressure time series from 30 June to 19 July 1975 
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Figure 9. RMS pressure errors and mean pressure errors for FNWC 
analysis at all AIDJEX stations for the period 1 May 1975 to 24 April 
1976. Dates indicated are AIDJEX days. Each point represents 20-day 
average value. The beginning of each month is indicated. RMS error 
is shown by a solid line, mean error by a dashed line. 
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Figure 10. RMS wind speed and direction difference between AIDJEX 
analyzed and FNWC analyzed and prognosticated winds. Lines marked 5 
differences for only those wind speeds greater than 5 m/sec. Points 
marked *G are aggregate differences for all four time periods. Time 
periods are A, 1 May to 20 May; B, 21 May to 9 June; C, 10 June to 
29 June; and D, 30 June to 19 July. 
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Figure 11a. Computed air stress magnitude (TAUA-MAG) from AIDJEX 
field (solid line) and FNWC field (dashed line) at AIDJEX main camp 
for days 15 to 25 May 1975. Units are dynes/cm2. 
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Figure lib. Computed air stress direction (TAUA-AZI) from AIDJEX 
field (solid line) and FNWC field (dashed line) at AIDJEX main camp 
for days 15 to 25 May 1975. Units are degrees true. 
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Figure 12a. 
field (solid 
15 to 25 May 



Computed air stress magnitude (TAUA-MAG) from AIDJEX 
line) and FNWC field (dashed line) at node 3,7 for days 
1975. Units are dynes/cm2. 
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Figure 12b. Computed air stress direction (TAUA-DIR) from AIDJEX 
field (solid line) and FNWC field (dashed line) at node 3,7 for days 
15 to 25 May 1975. Units are degrees true. 
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75RUN1F/P*=4E7/HBAR=300/G=0 /INITIAL STRESS P=.67 



Figure 13. AIDJEX computed air stress field at 1200GMT 18 May 1975. 

Scale vector represents 6 dynes/cm 2 . 
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FNWC RUN 135-2, USING REVISED FNWC AIR STRESS 



Figure 14. FNWC computed air stress field at 1200GNTT 18 May 1975. 
Scale vector represents 6 dynes/cm^. 
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X-POSITION VS Y-POSITION. CYCLE 0 TO 4981. NODE 7. 7 




Figure 15. Trajectories of AIDJEX main camp from 15 to 25 May 1975 
for AIDJEX simulation (solid line) and FNWC simulation (dashed line). 
( IS ) Starting position, ( # ) ending position, ( • ) actual measured 
position every six hours with AIDJEX dates indicated. Grid spacing 
is 5 km/div. 
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Figure 16a. X - Position time series of AIDJEX main camp as computed 
by AIDJEX simulation (solid line) and FNWC simulation (dashed line) . 
Ordinate spacing 5 km/div. 
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Figure 16b. Y - Position time series of AIDJEX main camp as computed 
by AIDJEX simulation (solid line) and FNWC simulation (dashed line) . 
Ordinate spacing 5 km/div. 
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Figure 17. Trajectories of node 3,7 from 15 to 25 May 1975 for AIDJEX 
simulation (solid line) and FNWC simulation (dashed line) . Grid 
spacing 5 km/div. 
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Figure 18a. X - Position time series of node 3,7 as computed by AIDJEX 
simulation (solid line) and FNWC simulation (dashed line) . Ordinate 
spacing is 5 km/div. 
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Figure 18b. Y - Position time series of node 3,7 as computed by AIDJEX 
simulation (solid line) and FNWC simulation (dashed line) . Ordinate 
spacing is 5 km/div. 
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75RUN1F/P*=4E7/HBAR=300/G=0 /INITIAL STRESS P=.67 



Figure 19. AIDJEX computed velocity field at 1200GMT 18 May 1975. 

Scale vector represents 30 cm/sec. 
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Figure 20. AIDJEX computed stress divergence field at 1200GMT 
18 May 1975. Scale vector represents 2 dynes/cm^. 
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FNWC RUN 135-2, USING REVISED FNWC AIR STRESS 



Figure 21. FNWC computed velocity field at 1200GMT 18 May 1975. 
Scale vector represents 30 cm/sec. 
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Figure 22. FNWC computed stress divergence field at 1200GMT 18 May 1975. 
Scale vector represents 2 dynes/cm^. 
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Figure 23a. Velocity magnitude of AIDJEX main camp as computed by 
AIDJEX simulation (solid line) and FNWC simulation (dashed line) . 
Speed units are cm/sec. 
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Figure 23b. Velocity direction of AIDJEX main camp as computed by 
AIDJEX simulation (solid line) and FNWC simulation (dashed line) . 
Direction units are degrees true. 
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Figure 24. AIDJEX computed air stress field at 1200GMT 22 May 1975. 
Scale vector represents 2 dynes/cm^. 
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Figure 25. AIDJEX computed velocity field at 1200GMT 22 May 1975. 
Scale vector represents 10 cm/sec. 
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Figure 26. FNWC computed air stress field at 1200GMT 22 May 1975. 
Scale vector represents 2 dynes/cm 2 . 
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Figure 27. FNWC computed velocity field at 1200GMT 22 May 1975. 
Scale vector represents 10 cm/sec. 
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